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Abstract 

We present a subtraction scheme for computing jet cross sections in electron- 
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QCD. In this first part we deal with the regularization of the doubly-real 
contribution to the NNLO correction. 
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1 Introduction 



QCD, the theory of strong interactions, is an important component of the Standard Model 
of elementary particle interactions. It is asymptotically free, which allows us to compute 
cross sections of elementary particle interactions at high energies as a perturbative expan- 
sion in the running strong coupling as{fiR). However, the running coupling as{fiR) remains 
rather large at energies relevant at recent and future colliders. In addition, to leading order 
in the perturbative expansion, the coupling varies sizeably with the choice of the (unphys- 
ical) renormalization scale fiR. In hadron-initiated processes, the situation is worsened by 
the dependence of the cross section on the (also unphysical) factorization scale fip, which 
separates the long-distance from the short-distance part of the strong interaction. Thus, 
a leading-order evaluation of the cross section yields rather unreliable predictions for most 
processes in the theory. To improve this situation, in the past 25 years the radiative cor- 
rections at the next-to-leading order (NLO) accuracy have been computed. These efforts 
have culminated, when process-independent methods were presented for computing QCD 
cross sections to NLO accuracy, namely the slicing [1112], subtraction O El E] and dipole 
subtraction [B] methods. In some cases, though, the NLO corrections were found to be 
disturbingly large, and/or the dependence on fi^ (and eventually ftp) was found to be stiU 
sizeable, thus casting doubts on the applicability of the perturbative predictions. When 
the NLO corrections are found to be of the same order as the leading-order prediction, the 
only way to assess the reliability of QCD perturbation theory is the computation of the 
next-to-next-to-leading order (NNLO) corrections. 

In recent years severe efforts have been made to compute the NNLO corrections to 
the parton distribution functions |7| and important basic processes, such as vector boson 
production |H1 IHl QHl ^| and Higgs production jHl El El E| in hadron collisions 
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and jet production in electron-positron annihilation These computations evaluate 

also the phase space integrals in d dimensions, thus, do not follow the process-independent 
methods used to compute the NLO corrections. 

The more traditional approach relies on defining approximate cross sections which match 
the singular behaviour of the QCD cross sections in all the relevant unresolved limits. 
Various attempts were made in this direction in Refs. [131201121112212312112312312111211 
l29] . In general, the definition of the approximate cross sections must rely on the singly- 
and doubly-unresolved limits of the QCD squared matrix elements. Although the infrared 
limits of QCD matrix elements have been extensively studied both at tree-level |H3 EH IH3 
E3I21I23E3EI1EHIE3I13I1I1, and at one-loop 1121 mm 113 Hi, the formulae presented 
in the literature do not lend themselves directly for devising the approximate cross sections 
for two reasons. The first problem is that the various single and double soft and/or coUinear 
limits overlap in a very complicated way and the infrared factorization formulae have to be 
written in such forms that these overlaps can be disentangled so that double subtraction is 
avoided. The second problem is that even if the factorization formulae are written such that 
double subtraction does not happen, the expressions cannot straightforwardly be used as 
subtraction formulae, because the momenta of the partons in the factorized matrix elements 
are unambiguously defined only in the strict soft and coUinear limits. In order to define 
the approximate cross sections one also has to factorize the phase space of the unresolved 
partons such that the singular factors can be integrated and the remaining expressions 
can be combined with the virtual correction leading to cross sections which are finite and 
integrable in four dimensions. 

In Ref. jlT] we presented a solution to the first problem, but did not explicitly define 
the approximate cross sections, which we left for later work. In this paper we turn to 
the second problem and define the complete approximate cross sections that regularize the 
doubly-real emission. For factorizing the phase space one has two options. On the one 
hand we may decompose the squared matrix elements into expressions that contain only 
single singular factors and use the factorization formulae as subtraction terms. In NLO 
computations this method was termed 'residuum subtraction'. On the other hand one may 
use exact phase-space factorization (keeping momentum conservation and particles on-shell) 
to maintain gauge invariance of the factorized matrix elements. In NLO computations the 
'dipole subtraction method' represents an example of this approach. 

The single singular factor decomposition cannot be followed in a NNLO computation 
because a singular factor in a doubly-unresolved region of the phase space naturally contains 
singular factors in singly-unresolved regions. Therefore, we have to use exact factorizations 
of the unresolved phase space measures. There are only two known ways of exact phase 
space factorization, one termed 'dipole factorization' j3, while the other called 'antennae 
factorization' Actually, both belong to the same general class. The important feature 
of these is that in order to factorize the unresolved phase space measures, two partons, 
called 'emitter' and 'spectator', are singled out for each subtraction term. The emitter 
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emits the unresolved parton and the spectator takes away the momentum recoil to maintain 
momentum conservation. In a NLO computation the choice for the emitter-spectator pair is 
naturally a colour-connected pair of partons in the soft factorization formula. The coUinear 
emissions are distributed among the soft ones using colour conservation. Although this 
trick provides a fairly elegant framework for computing NLO corrections, the unnatural 
distribution of collinear emissions seems impossible to maintain in a NNLO computation 
because it leads to simultaneously spin- and colour-correlated squared matrix elements for 
which collinear factorization formulae do not exist |T7].* 

In Ref. j2Hl the antennae factorization is used for NNLO subtractions, which is made 
possible by the use of colour-ordered subamplitudes, where the emitter-spectator 'antenna 
pairs' can naturally and unambiguously be selected because singular emission for a given 
amplitude can occur only between ordered pairs of momenta (the emitter and the specta- 
tor). Note however, that the effect of quantum interference in the squared matrix element 
mixes the colour subamplitudes in a rather complicated way and a general scheme to con- 
struct the subtraction terms in a process-independent way has not been given yet. It seems 
to us that in order to construct a general method for computing NNLO corrections we are 
forced to give up the antennae (or dipole) factorization of the phase space. 

In Ref. jlHl! we defined a new NLO subtraction scheme with exact phase-space fac- 
torization that can be generalized to any order in perturbation theory. In this paper we 
extend that scheme to computing the contribution of the doubly-real emission to the NNLO 
corrections. We demonstrate that the regularized cross section is indeed numerically inte- 
grable by computing the corresponding contribution to the cross section of electron-positron 
annihilation into three jets from the e'^e~ — >■ qqggg subprocess. 



2 Jet cross sections at NNLO accuracy 

The jet cross sections in perturbative QCD are represented by an expansion in the strong 
coupling as- At NNLO accuracy we keep the three lowest-order terms, 

^ ^LO I ^NLO I NNLO /O 1 ^ 

Assuming an m-jet quantity, the leading-order contribution is the integral of the fully 
differential Born cross section da^ of m final-state partons over the available m-parton 
phase space defined by the jet function J^, 

= [ da^Jm . (2.2) 



"Note a misprint in Ref. where 'soft' factorization is written in this respect instead of 'coUinear'. 
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The NLO contribution is a sum of two terms, the real and virtual corrections, 



J m+l Jm 



(2.3) 



Here the notation for the integrals indicates that the real correction involves m + l final- 
state partons, one of those being unresolved, while the virtual correction has m-parton 
kinematics. The NNLO correction is a sum of three contributions, the doubly-real, the 
one-loop singly-unresolved real-virtual and the two-loop doubly-virtual terms. 



Here the notation for the integrals indicates that the doubly-real corrections involve m + 2 
final-state partons, the real-virtual contribution involves m + l final-state partons and the 
doubly- virtual term is an integral over the phase space of m partons, and the phase spaces 
are restricted by the corresponding jet functions J„ that define the physical quantity. 

In d = A dimensions the two contributions in Eq. (|2.3|) as well as the three contributions 
in Eq. (j2.4j) are separately divergent, but their sum is finite for infrared-safe observables 
order by order in the expansion in as- The requirement of infrared-safety puts constraints 
on the analytic behaviour of the jet functions that were spelled out explicitly in Ref. [TF] . 

Following from kinematical reasons, fully inclusive observables can be accurately evalu- 
ated in QCD perturbation theory relatively simply. Since these observables are completely 
inclusive, no phase-space restriction has to be applied (J„ = 1 for any n). Real and virtual 
contributions can be combined at the integrand level resulting in the cancellation of soft 
and collinear singularities before performing the relevant phase-space integrations. Owing 
to these features, general techniques have been available for some time [52 E21 to carry out 
NNLO calculations in analytic form. 

QCD calculations beyond LO for inclusive quantities, such as jet cross sections or event- 
shape distributions, are much more involved. Owing to the complicated phase space for 
multiparton configurations, analytic calculations are impossible for most of the distribu- 
tions. Moreover, soft and collinear singularities are separately present in the real radiation 
correction (due to integrations over the phase space of the unresolved parton) and virtual 
contributions (due to integrations over the loop momentum) at the intermediate steps. 
These singularities have to be first regularized by analytic continuation in a number of 
space-time dimensions d = 4 — 2e different from four. This analytic continuation prevents 
a straightforward implementation of numerical integration techniques. 

The traditional approach to finding the finite corrections at NLO accuracy is to regu- 
larize the real radiation contribution by subtracting a suitably defined approximate cross 
section do"^'^ such that (i) do"^'^ matches the pointwise singular behaviour of da^ in the 
one-parton infrared regions of the phase space in any dimensions (ii) and it can be inte- 
grated over the one-parton phase space of the unresolved parton independently of the jet 




(2.4) 
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function, resulting in a Laurent expansion in e. After performing this integration, the ap- 
proximate cross section can be combined with the virtual correction do"^ before integration. 
We then write 



^NLO 



m+1 



+ 



d(j^ + / do-^+i 



Jm 1 (2-5) 



where both integrals on the right-hand side are finite in (i = 4 dimensions. The final result 
is that one is able to rewrite the two NLO contributions in Eq. ()2.3|1 as a sum of two finite 
integrals, 

= / da^+o + / da^^o , (2.6) 



J m+l Jm 

that are integrable in four dimensions using standard numerical techniques. 

The construction of the suitable approximate cross section da^'^ is made possible by the 
universal soft and collinear factorization properties of QCD matrix elements. Envisaging a 
similar construction for computing the NNLO correction, the universal infrared behaviour 
of the loop amplitudes and the infrared limits of the real-emission corrections at NNLO, as 
well as the singularity structure of the two-loop squared matrix elements has been computed 
1101 ED] • However, it is far more complex to disentangle these singularities at the NNLO 
accuracy than at NLO [Hj, thus up to now, process independent approximate cross sections 
for regularizing the da^^ and da^ terms have not been computed, but for the relatively 
simple case of e^e~ —>■ 2 and 3 jets, when the dependence on colour completely factorizes 
from all matrix elements 



In order to avoid this complexity in Ref. jool a new method has been developed for com- 
puting the QCD corrections by combining the real-emission and virtual corrections before 
integration for arbitrary jet function. The method is very simple conceptually. It considers 
the problem from a purely mathematical point of view: how to compute a complicated, 
but finite integral numerically? The first step is to map the phase spaces onto the unit 
hypercube of suitable dimensions. Then the singularities from inside the hypercube are 
removed to the edges of the cube by splitting appropriately the integrations and mapping 
them back to the [0, 1] interval. Next, the overlapping singularities are disentangled using 
sector decomposition j^H EHl EHl EZl EHl EH] • At this point the only factors in the integrand 
that lead to divergences are of the form A"^"*""^, therefore, the e poles can be extracted in 
terms of plus distributions The method is clearly completely general and its strength 
has already been demonstrated in various explicit computations [T^ll6|IT7| ll2j. Note how- 
ever, that the various mappings of the phase space as well as the sector decompositions are 
not unique. The particular choices depend on the analytic structure of the functions one 
has to integrate, namely, the squared matrix elements for the given process. In fact, the 
different terms in the squared matrix element may prefer different mappings as in the case 
of Ref. ^H] • This means that with this technique the construction of a universal program 
that requires only various matrix elements as input for computing NNLO corrections for 
arbitrary processes does not seem straightforward. 
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Such program exists for computing NLO corrections [HDlElinS] based upon the subtrac- 
tion method. Therefore, it is of interest whether the subtraction method can be extended 
to the computation of NNLO corrections. In the next section we rewrite Eq. (|2.4p such that 
each phase space integral is finite and thus can be performed numerically in four dimensions 
using standard Monte Carlo techniques. 



3 Subtraction scheme at NNLO accuracy 

Let us consider first the doubly- real contribution, dcx^^j- is divergent in the doubly- 
unresolved regions of phase space. In order to cancel the two-parton singularities we sub- 
tract the approximate cross section da^^'^^ that matches the pointwise singular behaviour 
of da^^2 dimensions in the two-parton infrared regions. Then we have 



a 



NNLO 



m+2 



+ / da 



1 RR,,A2 T 



m+1 



m+2 



+ 



J RR.Ao 



J yy 



(3.1) 



However the first integral is still divergent in the singly-unresolved regions of the phase 
space. In order to cancel these remaining singularities we subtract the approximate cross 
sections da^^';^^ and da^2^^ to obtain 



NNLO 



+ 



m+1 



j^RR 7 T /'Hrr^^'^i T 



m+2 '- 



da, 



RV 

m+1 



da. 



vv 



'm+2 



da 



RR,Ai 
m+2 



J, 



^'^m+2 



m+1 



1 RR.Aio 
'i^m+2 



1 RR,A^2 T 



(3.2) 



Here dcr^^'^^ and da^^'^^^ regularize the singly-unresolved limits of da^2 da^^'^^ 
respectively. For the construction to be consistent, we must also require that da^^'a ^ — 
do"^^'^^^ be integrable in the two-parton infrared regions of the phase space, ^ which restricts 
the possible forms of dcr^^'^^^ severely. In Eq. ()3.2|) the jet functions, multiplying each 
approximate cross section, organize the terms according to in which integral they should 
appear. In particular, da^^'^^ is multiplied by Jm+i, therefore, after integration over 
the phase space of the unresolved parton, it is combined with da^^^, while da^^2 is 
multiplied with J^, therefore, it is added back in the third line. The m + 2-parton integral 
above is now finite by construction. 

^'^m+2 1 which we have aheady taken into account in writing 



RR.A, 



T Formally this means that dcr^^j 
the expressions. 
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Next consider the real- virtual contribution da^_^i. It has two types of singularities: 
(i) explicit e poles in the loop amplitude and (ii) kinematical singularities in the singly- 
unresolved regions of the phase space. In Eq. (|3.2j) the former are already regularized. 
Indeed, unitarity guarantees that the second line of that equation is free of e poles if 
da^2^ is a true regulator of da^2 the one-parton infrared regions of the phase space, 
just as it does in NLO subtraction schemes. To regularize the kinematical singularities, we 

subtract the approximate cross sections da^^'f^ and ^/^ do"^^'^^ j'^^ which regularize the 

real- virtual cross section da^i and da^2^ , respectively, when a single parton becomes 
unresolved. Thus, the NNLO cross section is written as 



a 



NNLO 



m+2 



m+1 



1 PIR,A2 T 



da, 



RR,Ai 
m+2 



J. 



m+1 



da. 



RV.Ai 
m+1 



1 RR.A, \ Ai 
do-m+2 * 



e=0 
Jin. 



e=0 



+ 



d(T^+2 



1 R,R,.Ato 
dc^m+2 



+ 



da. 



RV,Ai 
m+1 



da. 



RR.Ai^Ai 
m+2 



(3.3) 



J rr 



£ = 



Since the first and second integrals on the right hand side of this equation are finite in 
d = A dimensions by construction, it follows from the Kinoshita-Lee-Nauenberg theorem 
that the combination of integrals in the last line is finite as well, provided the jet function 
defines an infrared-safe observable. 

The final result of these manipulations is that we rewrite Eq. ()2.4j) as 



a 



NNLO 



1 NNLO , 



m+2 



1 NNLO , 



m+1 



da 



NNLO 



(3.4) 



that is a sum of three integrals. 



1 NNLO 
dcrm+2 



1 NNLO 
'^'^m+1 



do'm+2'^^+2 



da^^ 

'^'^m+1 



1 FiR,,A2 J 



da 



RR,Ai 
m+2 



J m+1 



m+2 ^m,+l 
dc^m+1 + 



da 



RR.Ai^A 
m+2 



£=0 

J, 



£ = 



(3.5) 
(3.6) 



and 



dar+ 



1 RR.,A2 

*^^m+2 



1 R,R,,Ai9 

-dcT^+2 



da. 



RV,Ai 
m+1 



da. 



rr.aAAi 

m+2 



£ = 



Jm, (3.7) 



each integrable in four dimensions using standard numerical techniques. 



The subtraction scheme presented here differs somewhat from the one outlined in 
Ref. |17], where we assumed that da^^'^^^ can be defined such that (in the present no- 
tation) 



m+1 



1 RR,A2^2 
'^'^m+2 



da: 



RR,Ai^ Ai 
m+2 



J n 



finite 



(3.^ 
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in (i = 4 dimensions. However, as already emphasized in Ref. jUj, Eq. ()3.8j) does not follow 
from unitarity, rather it is a constraint on the definitions of da^^^ and da^2^^ . It turns 
out more convenient to drop this extra condition and rearrange the subtraction scheme as 
presented here. 

In this paper we present all formulae relevant for constructing dcrj^^g^ explicitly. The 
terms needed for defining dcr^^^° and dcr^'^'"^ will be given in separate papers. We use 
the colour- and spin-state notation introduced in Ref. [H]. The complete description of our 
notation can be found in Ref. |47j . 



4 Subtraction terms for doubly-real emission 

The cross section dcr^^g is the integral of the tree-level squared matrix element for m + 2 
parton production over the m + 2 parton phase space 

daZ2 = d<Pm-,2\M^Z2\\ (4.1) 
where the phase-space measure is defined as 

dUvi, . . . ,Pn; Q) = n J^Mv^^ (2^)'^^'^ [Q - ■ (4-2) 

In Eq. fl4.1|l (and all subsequent formulae) the superscript (0) refers to tree-level expres- 
sions. We disentangled the overlap structure of the singularities of |A1m+2p i^^o '^^^ pieces 

A2|A<!!!+2l', Ai|A<i||^2l' and Ai2|A<!^^2l' m Ref. gZI. These expressions are only defined 
in the strict soft and/or coUinear limits. To define true counterterms, they need to be 
extended over the full phase space. This extension requires a phase-space factorization 
that maintains momentum conservation exactly, but such that in addition it respects the 
delicate structure of cancellations among the various subtraction terms. 

The true (extended) counterterms may symbolically be written as 

da^^^^^ = d0^ \dv2\ J^^M'^Xf , (4-3) 
doT^t = d0™+i [dpi] Af^\M^Z2\' , (4.4) 

and 

daT^2'' = d<P^ [dp,] [dp,] 4?|-MSU|' . (4-5) 

where in Eqs. fl4.3|) - ()4.5p we used a formal, calligraphic notation (to be defined explicitly 

below) to indicate the extension of the terms A2|A1^+2p5 Ai|7\/l[J|l,„2p and Ai2|7W^+2p 
over the whole phase space that was written in exactly factorized forms, 

d0m+2 = dcpm [dp2] = dcpm+i [dpi] = dcpm [dpi] [dpi] (4.6) 
(the precise meaning of the factors [dpi] and [dp2] will be given below). 
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5 Singly- unresolved counterterms 



The singly-unresolved counterterm w4i|A^,^°L2p reads 



AiA^2U(M)r = E 



(5.1) 



Here all three terms are functions of the original m + 2 momenta that enter the matrix 
element on the left hand side of Eq. (j5.1|) . To shorten the notation we denote these momenta 
collectively as {p} = {pi, . . . ,pm+2}- Although the notation in Eq. ()5.1|) is very similar to 
the operator notation introduced in Ref. jlTj, it is important to understand that it is not 
meant in the operator sense, for instance, the last term on the right hand side does not 
refer to the coUinear limit of anything. Throughout this paper the subtraction terms are 
functions of the original momenta for which the notation inherits the operator structure of 
taking the various limits, but otherwise it has nothing to do with taking limits. 



5.1 CoUinear counterterm 



Counterterm 

The singly-collinear counterterm is 



1 



.2e /^j^{0) f[~\{ir) M £,(0) 



((0) ^r^T(«0 



where the kernels Pj^j^{zi^r, Zr,i, k±^i^r] £) are defined to coincide with the following specific 
forms of the Altarelli-Parisi splitting functions (valid in the CDR scheme) 

-1 + zf 



{r\P^l{z,,r,Zr,;e)\s)=5rsC, 



.2 



— ezr 



(5.3) 



{^i\P^%{z.,r,Zr,,k^^^,^,;e)W)=2CA 



kt „k 



2,r -L,2,r' 



\ 2 1 - e)zirZ., 



Zr,i Zi^r 



(5.4) 



.,i,r 



(5.5) 



In Eq. (|5.2j) the double superscipt on the left hand side means that on the right hand side 
of the equation both the matrix elements as well as the splitting kernels are at tree-level. 
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In Eq. (15 .Sj) we introduced our notation for the spin-averaged splitting function, 
The kernels are functions of the momentum fractions Zi^r and Zr^i that we define as 

~ - ^''^ and 7 ■ - ^'"'^ 7] 
dllU. I — , [O. i j 



y(ir)Q y{ir)Q 

where y{ir)Q = ViQ + VrQ with fjiq = 2pi ■ Q/Q'^, UrQ = 2pr ■ Q/Q"^ and Q^" is the total four- 
momentum of the incoming electron and positron. With this definition -|- Zr^i = 1. Note 
that the momentum fractions are nothing but the energy fractions of the daughter momenta 
of the splitting with respect to the energy of the parent parton in the center-of-momentum 
frame. The transverse momentum kj_^i^r is given by 

_ A M _ /■ M I /■. A. _ _ ^ _ _ Vir /p. o\ 

Ciiry{ir)Q Ciiry{ir)Q 

Here yir = 2pi ■ Pr/Q'^ while pf^ and a^r are defined below in Eqs. (jS.llj) and ()5.12j) respec- 
tively. This choice for the transverse momentum is exactly perpendicular to the parent 
momentum p'^^ and ensures that in the coUinear limit Pi\\Pr, the square of fc^ ^ ^ behaves as 

as required (independently of Qr)- In a NLO computation this feature is sufficient to 
ensure the correct collinear behaviour of the subtraction term. In a NNLO computation 
in addition to Eq. ()5.9|) it is also important that k'^^^ itself vanishes in the collinear limit, 
and it is that convenient that k'^^^ is perpendicular to Q^^. These conditions are fulfilled if 
we choose ^ 

Cir ~ (yZr,i Zi,r) ■ (5.10) 

'^iry'i^Q 

With this choice k'^^^ ^ k'^^ in the collinear limit as can be shown by substituting the Su- 
dakov parametrization of the momenta into Eq. (j5.8|) (with properly chosen gauge vector). 

Momentum mapping and phase space factorization 

The m + 1 momenta, {p}m+i = {Pi^ ■ ■ ■ ^Pir^ ■ ■ ■ ,Pm+2}-, entering the matrix elements on 
the right hand side of Eq. ()5.2|) are defined as follows 

Pt = Y^(P^+Pr-»irQn, P'n=T^P'n, n=^t,r, (5.11) 



where 

_ 1 r 



y^Q- ^yf^r)Q-^y^r \ ■ (5-12) 

^li k'^ ^ J, does not vanish in the coUinear hmit then the iterated coUinear-triple coUinear counterterms 
of Sect. [Tn] do not have the correct (strongly-ordered) coUinear behaviour. 
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The total four-momentum is clearly conserved, 

m m 
n n 

For further convenience let us denote the momentum mapping introduced above as 

M ^ {piSi ' (5-14) 

where i and r are any two labels of momenta that appear on the left hand side, but not on 
the right hand side. 

The momentum mapping of Eq. (jS.lH) leads to exact phase space factorization in the 
form 

; Q) [drf5+i(Pr,Pir; Q)] , (5.15) 

where, as indicated, the m+1 momenta in the first factor on the right hand side of Eq. ()5.15|) 
are exactly those defined in Eq. ()5.11|) . The explicit expression for [<^Pi^m^i{Pr,Pir] Q)] reads 

[drf5+l(Pr,Pir-; Q)] = jS+l{Pr,Pir; Q) ^^^^^^l^ > (5-16) 

where the Jacobian is 

r^ fl - a- Bfl - a- ) 

jS^liPr,P^r; Q) = VTrQ 7 ^ I ''^ ■ (5.17) 

^{^-yTrQ)(^ir + yr^ +yTrQ-yrQ 

In this equation a^r is the physical (falling between and 1) solution of the constraint 



In order to implement the subtraction scheme, this solution is not required, the momentum 
mapping given in Eq. (jS.llj) is sufficient. 

The collinear momentum mapping of Eq. ()5.11|) and the implied phase-space factor- 
ization of Eqs. ()5.15p - ()5.17|) are represented graphically in Fig.Q The leftmost picture 
represents the (m + 2)-parton phase space d0m+2({p}; Q)- In the circle we denote the num- 
ber of final-state partons. The picture in the middle represents the result of the mapping 
of momenta in Eq. (jS.llj) . The dots between the momentum pi^ and the circle with the two 
momenta Pi and Pr means that the latter two are replaced with pi^. This mapping implies 
the exact factorization of the phase space, written in Eq. ()5.15|) and represented by the pic- 
ture on the right. The first factor is the phase space of (m + l)-partons, d0m+i({p}^+i; Q), 
and the second is [dp^*^'']. In the latter the box represents the Jacobian on Eq. ()5.17p and 
the line means the one-particle phase-space measure. We shall use similar graphical rep- 
resentations of other momentum mappings and implied factorizations of the phase space. 
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Figure 1: Graphical representation of the singly-colhnear momentum mapping and the 
imphed factorization of the phase space. 

In the colhnear hmit, when pfHp^' ^^e transverse momentum behaves as in Eq. ()5.9|) . 
Zi,r Zi and Zr^i —>■ z^- Furthermore, tends to zero so f)^^ —* Vi + Vr pn Pn, 
i.e. the tildes disappear from Fig.d Consequently, the counterterm reproduces the colhnear 
behaviour of the squared matrix element in this limit. 

5.2 Soft-type counterterms 

Counterterms 

The soft-type^ terms are the singly-soft and singly soft-coUinear counterterms 

Si'fHip}) = -8vras^^^$^E^'5^^-WI-^!^li,{..)({^}^ii)l'' (5-19) 

C.rSi:^'\{p})=8na.l^'^-^T^.\M^ZA^^^^^ (5.20) 

Sir Zr,i 

If r is a quark or antiquark, Sr^''^\{p}) and Cj^iS^°'°^ ({p}) are both zero. The eikonal factor 
in Eq. (ICTHl is 

5..(r) = ^, (5.21) 

and the momentum fractions entering Eq. (|5.2Up are given in Eq. (j5.7|) . The operator Tj is 
the colour charge of parton i jH]. 

Momentum mapping and phase space factorization 

The m + 1 momenta, {pj^l^^ = {pi, . . . ,Pr-i,Pr+i, ■ ■ ■ ,Pm+2} (the momentum with 
index r is absent), entering the matrix elements on the right hand sides of Eqs. fl5.19|) 
and (j5.2(Jj) are defined by first rescaling all the hard momenta by a factor 1/A,. and then 

§Thc expression 'soft-type' refers to the momentum mapping used to define these terms. 



13 



transforming all of the rescaled momenta as 

= A^[g, (g - Pr)/KWJK) , n^r, (5.22) 



where 
and 



K = ^/i^y^. (5.23) 

AnA-,X|^.^- ^'^--™:^W gj^. (5.24, 

The matrix A^[i^, K] generates a (proper) Lorentz transformation, provided = K ^ 
0. Since is massless (p^ = 0), the total four- momentum is again conserved, 

m+l m+1 

Q'=P^.+Y.Pn = Y.Pn- (5.25) 

n n 

We will find it convenient to introduce the notation 

M ^ {P}t\i (5.26) 

to denote the above momentum mapping. Here r is the label of any momentum of the 
original momentum set {p} that is absent from the set on the right hand side. 

The momentum mapping of Eq. ()5.22j) also leads to exact phase space factorization 

d^m+2i{p};Q) = d^m+imtU Q) Wl-m+M-; Q)] ■ (5.27) 

The m + l momenta in the first factor on the right hand side are those of Eq. fl5.22|) . The 
factorized one-parton phase space [dp[r^_,_]^(pr.; Q)] is 

[dMi-M(p.; Q)] = Jtl^M; Q) j^Mp") > (5.28) 

with Jacobian 

JS^^l{Pr■,Q)=K^'-'^-'Qi^r). (5.29) 



Similarly to the graphical representation of the coUinear momentum mapping and 
phase-space factorization, depicted in Fig.^ we present the graphical representation of 
Eqs. fj5.22|) - ()5.29|) in Fig.|21 The middle picture in this figure represents the soft-type map- 
ping of Eq. ()5.22j) . The dots between the main circle and the circle with momentum pj. 
mean that the latter does not take away momentum from (see Eq. fl5.25|) ). 

The soft-type terms are defined on the same phase space, therefore, in the coUinear limit, 
when Pi\\p!^, the soft-collinear counterterm regularizes the kinematical singularity of the 
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Figure 2: Graphical representation of the singly-soft momentum mapping and the implied 
factorization of the phase space. 

soft counterterm by construction. In the soft limit, when — > 0, in Eqs. ()5.1H) and ()5.12|) 
air —>■ and pf^ — > pf, therefore, the momenta obtained in the collinear mapping tend to 
the same momenta as those obtained in the soft mapping, i.e. all momenta with tilde tend 
to the corresponding original momenta both in the case of collinear and soft momentum 
mappings, with the soft momentum being dropped as shown in Fig. El Thus the soft- 
coUinear counterterm regularizes the kinematical singularity of the collinear counterterm. 
At the same time the soft one regularizes the squared matrix element by construction. 



i i 1 




Figure 3: Graphical representation of the soft limit of collinear and soft-type mappings. 

In closing this section, we note that the jet function Jm+i that multiplies da^^'^^ in 
Eq. (|3.3|) is a function of the m + l momenta {p}^^li and for the collinear and 

soft-type countertems, respectively. 
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6 Doubly-unresolved counterterms 



The doubly-unresolved counterterm is 



6 



+ E IC,^;A'''\{P})\ +l^l's'\{p})] ■ (6.1) 



6.1 Triple collinear counterterm 

Counterterm 

The triple collinear counterterm reads 

Cfrf\{P}) = i8^C^sf^'n'^{M'^\{p}t''^)\PMrf^^^^^^ ■ 

(6.2) 

The Pfif^f^{{zj^kh Sjk, k± j^i:i}; e) kernels are defined to be the specific forms of the triple 
parton splitting functions introduced in Ref. jl7j, with one important modification: In the 
gluon splitting functions Pg^q^q^ and -Pgig^g^', the azimuth-dependent terms that depend on 
the transverse momenta have always to be written in the form k'^jk'^ k/^^,i ' ^-L,^ t)e 
equal to j), otherwise the collinear behaviour of the counterterm cannot be matched with 
that of the single collinear counterterm in the singly-unresolved phase space region. The 
correct azimuth-dependence can be achieved by making use of the following identites: 



k± j^l j = ( ~ ~ Zj)SjM + ZjSkl ) — - , (6.3) 

/ \ k^ k^ 

2 K,jkl,k = [^jk + 2 ZjZkSjki - ZiSji^ik) - ZjSi(^jk)J ^ ^'^ ^'^ , (6.4) 

where {A;, /} = {i, r, s} \ {j} and j can be i, r or s. Note that Ref. does not consider the 
splitting function for the case of final-state fermions with identical fiavours. The reason is 

^In the case of the g ggg sphtting, Ref. 07] uses the original definition of Ref. j^T], not spelled out 
explicitly. 
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that it consists of two terms of the type of different flavours plus a third term correspond- 
ing to an interference contribution denoted by (Pq^lq^) in Ref. jHZj. The different flavour 
contributions were given in Ref. jl7j, while the interference term can directly be taken from 
Ref. with the simple substitution of indices 1 ^ r, 2 ^ s and 3 ^ z in order to match 
our notation. This latter term does not require any special care because it does not have 
a leading singularity in any of the singly-, or other doubly-unresolved regions of the phase 
space, apart from the triple collinear one. 

The momentum fractions in Eq. ()(i.2|l are defined similarly as for the collinear case given 
in Eq. (jKTjl 



y{irs)Q y{irs)Q y{irs)Q 

with Zi^rs + Zr^is + ^s,ir = 1- The transverse momentum k±^r,is is 

^1,r,is ~ Cr,isPi ~ Ci,rsPr + Cr,isPs ~ Cs,irPr + CrisPirs ) (6.6) 

where" 

> y{rs)i ^ Vir Vrs '^^r^isVirs in y\ 

^irs y(irs)Q ^irs y^j^Q 

The expression for k^^g^i^ is obtained from Eq. (j6.6p by simply interchanging the indices r 
and s, while k^^i ^s = —k±^r,is — k±^s,ir- Similarly, (r,is and (s,ir are obtained from Q^rsj while 
(irs and (sir from (ris by interchanging the indices. We define pf^^ and airs in Eqs. fl6.9|) 
and ()6.10|) . This choice for the transverse momentum is also perpendicular to the parent 
momentum p^^^ as the one given in Eq. (jS.Sj) and it ensures also that in the triple collinear 
limit pf relations of the type (see Ref. [SZj) 

Srs — Zr,is Zs,ir I I ; v'-''"/ 



as well as k±^r,is k±^r are fulfilled. In a NNLO computation the longitudinal component 
in Eq. ()6.6p does not give any contribution due to gauge invariance, therefore, we may 
choose Cris = 0. 

Momentum mapping and phase space factorization 



The matrix elements on the right hand side of Eq. ()6.2j) are evaluated with the m 
momenta {pjm''^ = {pi, . . .,Pirs, ■ ■ ■,Pm+2} defined as 

pL = y^ ip': + P'f+P':-c^^rsQn, Pn = T^ P^ n^t,r,s, (6.9) 



with 

_ 1 



y{irs)Q - sJylrs)Q-^yirs ■ (6.10) 



"Note that the indices of Cris are ordered! 
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Figure 4: Graphical representation of the triply-coUinear momentum mapping and the 
imphed factorization of the phase space. 

Clearly Eqs. (j6.9p and (|6.1Up are generalizations of Eqs. (|5.11|) and (j5.12j) . The total four- 
momentum is again conserved, 



The momentum mapping of Eq. ()6.9p leads to exact phase-space factorization in a form 
very similar to Eq. (jS.lSj) 

d0^+2({p}; Q) = <lM{pYr^; Q) Wt^\Vr,Ps,P^rs■, Q)] , (6.12) 

where the m momenta in the first factor on the right hand side of Eq. (j6.12j) are those given 
in Eq. 1)6. The explicit expression for [(^p^^m {Pr-,Ps-,Pirs'-,Q)] is 

[^ptm{Pr,Ps,Pirs\ Q)] = J^^\Pr,Ps,Pirs] Q) (^~^'^+*^^^) ' (^'^^^ 

where the Jacobian is 

j""Hv P P ■Q)-v~ (1 - 9(1 - a,..) 

^l-"- y irsQl^irs ^ y(r-s)irs ^ y irsQ y{rs)Q 

In this equation airs is the physical (falling between and 1) solution of the constraint 
P^ 

^ = (1 - VI^sq) «L + (2/(r.)1^ + Vlir.Q - yirs)Q) O^irs + Vrs " = . (6.15) 

We present the graphical representation of Eqs. ()6.9p - (j6.14j) in Fig.0] The meaning of the 
various graphical elements is analogous to those in Fig.[TJ 

In the triply-coUinear limit, when pf||Prlbs5 ^^e transverse momenta behave as in 
Eq. (j6.8|) . Zi^rs Zi, Zr^is — Zj. and Zg^ir — > Zg. Furthermore, airs tends to zero so 
Pirs ~^ Pi + Pr ~'<~ Ps Pn ^ Pn (the tildcs disappear from Fig.|H). Consequently, the 
counterterm reproduces the coUinear behaviour of the squared matrix element in this limit. 
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6.2 Double collinear counterterm 



Counterterm 

The double collinear subtraction term reads 

= (8vras/i'^)2^ (6.16) 

where Pi^\{zk, zi, k±^; e) are given in Eqs. (j5.3p - (j5.5|) . The momentum fractions are defined 
in Eq. (jo. 7j) . and the transverse momentum k±^ir-js is 

kl,r■Js = C^,rP^-Cr,^p': + C^rPfr, (6-17) 

i.e., it is formally identical to the k'^ ^ ^ defined in Eq. ()5.8j) . although the definition of pf^ is 
different in the two cases (cf. Eqs. ()5.11|) and ()6.18p ). In the NNLO computation we may 
choose (ir = in Eq. ()6.17|) . Of course, k±js;ir given by Eq. ()6.17|) after the interchange of 
indices {i ^ j,r ^ s). Eq. (j5.12j) defines a^r and ajg, the latter with the same change of 
indices as before. 

Momentum mapping and phase space factorization 

The m momenta {p}m'''^^ = {p, ■ ■ ■ -.Pir-, ■ ■ ■ ,Pjs-, ■ ■ ■ ■,Pm+2} entering the matrix element 
on the right hand side of Eq. ()6.16p are again given by a simple generalization of Eq. ()5.11|) 

P'r = +Pr- ' P's = T—-^-—-iP' +Ps- ^,sQn , 

Pn = 7-z — -Z Pn^ n^i,r,j,s. (6.18) 

The total four-momentum is clearly conserved. 

The momentum mapping of Eq. ()6.18|) leads to the following exact factorization of the 
phase space: 

d<Pm+2{{phQ) = d<j>U{p}t'''^'^; Q) [dp^^'\pr,Ps.P^r,PJs■, Q)] • (6.19) 

where the m momenta in the first factor on the right hand side of Eq. ()6.19|1 are given by 
Eq. dnUI) and [dp^''^^'^^ (p, 

iPsyPiryPjs'iQ)] reads 

; Psy Piry Pjsj Q)] = J^T'iPr ! Ps^ Piry Pjsj Q) 

x^MpD^MpD- (6-20) 
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Figure 5: Graphical representation of the doubly-colhnear momentum mapping and the 
imphed factorization of the phase space. 



In this equation the Jacobian factor can be written as, for instance, 



+ 
+ 



a 



+iyi^Q - yrQ){yrsQ - y^Q) - (yi^Q + yrQ)ysrs - (yrso + y^Qh 



'jsQ 



SJS 



(6.21) 



where air and ajs are the physical (both are positive and their sum falling between and 
1) solutions of the coupled constraints 



p. 

Q 



P 

^ = (1 - yjsQ) "jl - VTsQ^i^^l^ + ^ysjs + yjsQ - ysO) + VsTs^is - 1/, = . (6.22) 



SJS 



The analytical solution of these equations is rather complicated. In a numerical calculation 
one can always find the solution numerically. However, for implementing the subtraction 
scheme, we need only the momentum mapping given in Eq. ()6.18p . In order to integrate the 
subtraction term over the factorized phase space [dp2^^^\pr^Ps,Pir^Pjs] Q)]^ we can choose 
air and ajs as integration variables and thus the solution of the coupled quadratic equations 
can be avoided. We present the graphical representation of Eqs. ()6. 18^ - 1)6. 211) in Fig. El 

In the doubly-collinear limit, when and p^j\\p^-, the transverse momenta behave 

as in Eq. Zk,i — ^ Zk {k, I = i, j, r, s). Both air and ajs tend to zero so pf^ Pi +Pry 
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Pjs ~* Pj + Ps-: Pn ~^ Pn- Consequently, the counterterm properly regularizes the squared 
matrix element in this limit. 



6.3 Double soft-collinear-type counterterms 

Counterterms 

We refer to the three terms on the second line of Eq. ()6.H) as double soft-collinear- 
type terms because they are all defined using the same momentum mapping, that of the 
soft-coUinear subtraction term. Explicitly these terms read 

(M) = -(8^«s/x^^)'EE^'5^'^(^) 

^^{M^:H{p}li'''^)\T,T,Pf;i{z,^^^ (6.23) 

°{ir)s ^s,ir 

^-{M^:H{p}li''^)\P^%i^^,r,Z^^^^^^^ (6.24) 

(M) = {S^a^f^'n'^'f^T^, 

><-{M^^\{p}ii'''^^)\Piliz,,r,Zr^^^^^^^ (6.25) 

Sir 

Here we remind the reader of what we wrote below Eq. ()5.ip about the notation. For 
instance, on the left hand side of Eq. ()6.24|) Q^^GSj^.f^ denotes a function that is defined 
by the function of the original momenta given on the right hand side. The eikonal factor 
entering the soft-coUinear counterterm GSjY.f'*({p}) is given in Eq. 1)5.211) . Nevertheless, we 
record explicitly that whenever j or / in Eq. (|6.23|) is equal to (ir), Eq. ()5.2ip evaluates to 

S(ir)i{s) = = — . (6.26) 

The momentum fractions appearing in Eqs. ()6. 23^ - 1)6. 25j) are the same as those for the 
singly-coUinear counterterm (see Eq. ()5.7)l ). while transverse components k±^i^r are the im- 
age of the transverse momentum /c_L,i,r in Eq. ()5.8|) under the soft map of Eq. ()5.26p 

^±,i,r ^ ^±,2,r • (6.27) 

Explicitly we have 

~kl„r = KiQ^ (Q - P^)/M{kl,,/\,) , (6.28) 
21 



where and the matrix are defined in Eqs. ()5.23jl and ()5.24j) . The hatted momentum 
ps is given in Eq. (|6.3(jp with n = s. 

Momentum mapping and phase space factorization 

The m momenta {p}m^^^ = {pi, ■ ■ ■ ,Piri ■ ■ ■ ^Pm+2] {Ps is absent) entering the matrix 
elements on the right hand sides of Eqs. ()6.23|) - ()6.25|) are defined by the composition of a 
single coUinear and a single soft mapping as follows: 

{P} ^ {P}^i . (6.29) 

The first mapping is a coUinear-type one leading to the hatted momenta 

p^ = Y^{pf + p'^-(^^rQn, p;^ = -^— pt:, n^i,r. (e.ao) 

where air is given in Eq. fj5.12j) . followed by a soft mapping 

= A^AQ, {Q - ps)/Xsm/Xs) , n^s, (6.31) 

where is given in Eq. ()5.23p (with r — > s). The other of the soft and coUinear mappings 
is not crucial. The order chosen here leads to simpler integrals of the singular factors over 
the unresolved phase space. 

This momentum mapping leads to exact phase-space factorization 

d<Pm+2i{p};Q) = d^mmii^''^; Q) Wtl+l{Pr.P^r■. Q)] Wfm{Ps; Q)] • (6.32) 

The coUinear and soft one-particle factorized phase space measures [dpi^^j^i{pr,Pir]Q)\ 
and [(lpi}i^{ps]Q)\ are given respectively in Eqs. (|5.16p and (j5.28|) . We show a graphical 
representation for Eqs. (|6.29p - (|6.32p in Fig. El 

The counterterms in Eqs. ()6.23|) - ()6.25|) are defined on the same phase spaces. There- 
fore, in the triply-coUinear limit, when the term in Eq. ()6.24|) regularizes that 
in Eq. ()6.23p . while Eq. ()6.25|) is integrable. Conversely, in the doubly-coUinear limit, when 
p^WPr and p^Wp^, the term in Eq. (j6.25p regularizes that in Eq. (j6.23p . while Eq. ()6.24p is 
integrable. In the soft-coUinear limit, when PiWp^ and — > simultaneously, both the 
triply-coUinear mapping, defined in Eqs. (|6.9p and (|6.1(J|) . and the doubly-coUinear map- 
ping, defined in Eq. ()6.18p . approach the corresponding limit of the iterated mappings of 
Eqs. fl6.29|) - ()6.31|) . This limit of the momentum mappings is depicted graphically in Fig.[7| 
We then conclude that the subtraction term in Eq. ()6.24j) regularizes the kinematical sin- 
gularities in Eq. (j6.2j) . while Eq. (j6.25p regularizes the singularities in the doubly-coUinear 
subtraction in Eq. In this limit, the subtraction term in Eq. (j6.23p is a local coun- 
terterm to the squared matrix element by construction. 



22 




m+2 





m+2 











{ "1 ) ir ® 


("■) 


— r (g) 





m+2 



Figure 6: Graphical representation of the soft-colhnear momentum mapping and the im- 
phed factorization of the phase space. 
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Figure 7: Graphical representation of the soft-coUinear limit of the triply-, doubly- and 
soft-collinear momentum mappings. 
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6.4 Double soft-type counterterms 



All terms on the last two lines of Eq. (jfi.lll are constructed using the same momentum 
mapping as the doubly soft counterterm Si^s'^\{p}) itself. We refer to these as doubly 
soft-type terms. 

Counterterms 

We begin by exhibiting the soft-type terms explicitly for the case of double soft gluon 
emission and then present counterterms to it in the various doubly-unresolved regions. We 
have 



'5S(M) = (8vra,/i2^)^ 



i,j,k,l 



i.k 



(6.33) 



A 













1 



+ 



1 



(1 £) (^Si^Zg^,/. SigZf^ig^ 
2^2 



[r ^ s] 



(6.34) 



Sir 



c cs >s(°'°) 

ITgSg ^fg^Sg TgSg 



2e\2 



TiTi\Mt\{PY::^)? ■ 



(6.36) 
(6.37) 



The momentum fractions were defined in Eqs. ()5.7|1 and (j6.5|l . The two-gluon soft function 
appearing in Eq. ()6.33|1 is 



I ^ir^ks H" ^is^kr 
Si(rs)^k{rs) 



—S.kirs) , (6.38) 
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where 

^^"•^(r, s) = S,kis) (SUr) + S^r) - S^r)) (6.39) 

is the form of this function in the strongly-ordered approximation and Sik{rs) is given by 
Eq. dOU, 

S^rs) = ■ (6.40) 

The discussion below Eq. fl6.25|) about the eikonal factor appering in the double soft- 
collinear subtraction term (Eq. (j6.23p ) and especially Eq. fj6.26p apply also to the soft- 
coUinear limit of the doubly-soft subtraction, Eq. ()6.35|) . 

The only nonzero terms for the case of emission of a soft quark-antiquark pair are 
Sr^'sg{{p}) and C^rgSg^r^'^q {{p}) ^ the remaining terms all vanish. Explicit expressions for 
these nonzero terms are 

'5i°J(W) = (8vra./i^^)^-^TR 

^rs 



X 



E E I ^ - 2^ \M^Z.^^mt'^)\' , (6.41) 



Si{rs)Sk(rs) ^i(rs) 



\M^^\{p}t'^)\' ■ (6.42) 



Momentum mapping and phase space factorization 

The m momenta {p}m^^ = {pi, ■ ■ ■ ,Pm+2} {Pr and Ps are omitted from the set) entering 
the matrix elements on the right hand sides of Eqs. (j6.33p - ()6.37p and Eqs. ()6.4ip - ()6.42p are 
given by a generalization of the singly-soft momentum mapping of Eq. (j5.22j) to the case of 
two soft momenta, 

= A^AQ, {Q-Pr- Ps)/A..](K/Ar.s) , n ^ r, s , (6.43) 

where 

Ks = \l 1 - (y(rs)Q - Vrs) (6.44) 

and A(^[(5, {Q — Pr — Ps)/'^rs\ is the matrix given in Eq. ()5.24|) . This momentum mapping 
again conserves total four- momentum. 

The m + 2 parton phase space factorizes exactly under the mapping of Eq. ()6.43|) . We 
find 

d0™+2({p}; g) = d0^({p}(;^); g)[dpg(p,,,p,; Q)] , (6.45) 
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Figure 8: Graphical representation of the doubly-soft momentum mapping and the implied 
factorization of the phase space. 

where the m momenta in the first factor on the right hand side are those defined in 
Eq. (j6.43|) . The explicit form of the factorized two-parton phase space is 

[dpS(p.,P.; Q)] = jSiPr,Ps; Q) ^^Mv") -(^MvD ' (6-46) 

with Jacobian 

jS^Pr.Ps: Q) = XT}'-'^-'eiXrs) . (6.47) 

We show a graphical representation for Eqs. ()6.43|) - ()6.47|) in Fig.|H| 

The counterterms in Eqs. ()6.33p - ()6.37|) are all defined on the same phase spaces. There- 
fore, all the cancellations among these terms in the triply-collinear limit, when Hp^HPs' 
in the doubly-coUinear limit, when pf\\pif and PjWPg, and in the soft-coUinear limit, when 
PiWPr and — 0, take place in just the same way as for the QCD factorization formulae 
described in Ref. jl?]. Consequently, the combination of the terms in Eqs. ()6.33p - ()6.37|) as 
present in Eq. ()6.1|) is integrable in four dimensions in all of these limits. The same is true 
for the difference of the two terms in Eqs. (|6.4ip and (|6.42p . 

In the doubly-soft regions of the phase space, when and 0, the momentum 

mappings in Eqns. (j6.9|) . (j6.18j) and in (j6.29p approach the corresponding limit of the double 
soft-type mapping in Eq. ()6.43|) . shown graphically in Fig. El Therefore, the cancellation of 
kinematical singularities among the various subtraction terms in Eq. ()6.H) takes place in 
just the same way as for the QCD factorization formulae described in Ref. jlTj, with the 
exception of the double soft subtraction terms in Eqs. (j6.33|) and (j6.41|) . In this limit, the 
latter provide local counterterms to the squared matrix element by construction. 
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Figure 9: Graphical representation of the doubly-soft limit of the triply-, doubly-, soft- 
collinear and doubly-soft momentum mappings. 



7 Iterated singly-unresolved counterterms 



In Ref. IHj, we introduced the term Ai2|A1^1^2p that reproduced simultaneously both the 
singly-unresolved limits of the doubly-unresolved limits and the doubly-unresolved limits 
of the singly- unresolved limits of the squared matrix element. The structure of these terms 
was such that Ai2|A^^)^2p '^^^ defined to be A1A2I A^^!^2p- When extending the terms 
Ai|7Vl^|,_2|^ and A2|7Vl^+2l^ o^^r the whole phase space to obtain the subtraction terms 



Ai|Al||^^2p ^^cl »4.2|7Vi^+2p5 ^6 defined the momentum mappings such that this structure 



,(0) 



could be preserved, i.e., the iterated singly-unresolved counterterm Ai2|A1m+2p is just the 



,(0) 



singly-unresolved subtraction for the doubly-unresolved counterterm A2|A1^+2l 
formally we can write 



Thus 



Ai2\M^Z2m)\' = Yl 



(7.1) 
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where the three terms in Eq. ()7.H) each evaluate further into the expressions,** 



r^k,t 



^kt^ktr ' ^kt^kt;r ^ kt^ ktr^ kt;r ^kt^rkf^kt 



,(0,0) 



(0,0) 



,(0,0) 



, (\r, ^(0,0) _r, p nAOfl) \ 

2-^ I 2 ^kt^ir\kt ^kt^ir;kt kt;r I 
i^r,k,t \ / 



I r o(o,o) 

"T ^kt'^kt ' 



(72) 



,(0,0) 



(0,0) 



o r o(o,o) <^ <^(o,o) , <^ ^ <,(o,o) 



5,5,?'°) ^ (7.3) 



and 



r^k,t 



'^kt^t'^krt 



/7 C /7 <,(0,0) _ ^ <^ ^ <^(0,0) , <-(0,0) 
''^kt'^t^kTt'^rt ^kt'^&rkt'^kt ^kt'^t'^rt 



^kt t rkt kt 



+ C,,S,Sf^'\ (7.4) 



In this section we spell out explicitly all subtraction terms that appear in Eqs. ()7.2p - (j7.4p . 
These subtraction terms are defined on factorized phase spaces that are obtained using two 
singly- unresolved coUinear and/or soft-type mappings iteratively. 



7.1 Iterated collinear counterterms 



7.1.1 Iterated collinear— triple collinear counterterm 



Counterterm 



The counterterm corresponding to the collinear limit of the triple collinear subtraction 



IS 



CktCktr\{v}) = {8Tia,fi 



2e\2 



1 1 



Skt S 



(7.5) 



kt f 



**At the level of the factorization formulae in Ref. 07] we neglected subleading terms in the triply- 
collinear limit of the soft-coUinear formula (see Eqns. (4.33) and (4.34) in 07|), which we do not apply 
here. Therefore, the soft limit of the doubly-real subtraction is slightly different from that in Ref. |47| . 
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where j, I = k,t and j, I = kt , f, defined in Eqs. ()7.12|) - ()7.14j) . Note that in our convention 
the ordering of the labels on the splitting-kernels is usually meaningless, but in the strongly 
ordered kernels P^f ^jjj} the ordering matters. As a result, the same triple-parton splitting 
function may have different strongly-ordered limits, which can be distinguished by the 
momentum labels in the kernel, once the ordering of the limits is fixed by the momentum 
mapping, given in Eq. ()7.12|) . Thus, in the following kernels, we do not further mark 
the abelian and non-abelian limits as we did at the level of the factorization formulae in 
Ref. For quark splitting we have 

5S.0. (0) 



' f,kt ' L,kt ,f' 



T, 



R 



k,t^ kti 



(^1 A'sfH^^fc,*, zt,k, z^^., 2;,^^ ; e)\s') = 



(7.6) 



(7.7) 



^Qfl^^ ^^r, kt ' ^kt ,f ' ^) 



X I — — 1 + Cf(1 - e)Zk,tZt,k I z 



'^L,k,t^ ktr 



Zt,k Zk^t 

where Sfkj_ k t — '^Vf ■ k±^k,t- For gluon splitting we define 

i^l^grqkqt i^k,t, Zt,k, k±^k,t, z-g.^^, z~^-g_ , k^ -g^y, e)\v) = 



{7.t 



(7.9) 



2Ca7r 



9 I \- + Zk,tZt,k ,^2 



^ kt ^f,kt 



'^±,k,t^ktr 



4:ZutZt,k 



■^kt,f '^±,k,t'^±,k,t 

^f,M ^^,k,t 



-ACa{1 - e)z.^'^^z-g^^.P^fljXzk,t,zt^k,k±^k,t;e) 



A-,f, kt 



(7.10) 
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and 



^r, fct ^ ^ kt ,f 



2^fc,t , Zt^k\ ^ ay 1 — e ^fk^^k,t o/l \ ^Tt ,r k±,k,tkl,k,t 
— + —]+9^ ^k,tZt,k 2 - 2(1 - £)^A:,i2:t,fc— -— -2 



,k,f^ktf '^f,kt '^±,k,t 

K'ktK-kt 

-ACA{l-e)z. jrtZTPt^,P^'ll{zk,uZt,k,k^,k,t,£) ^'p' j^''"'"'* • (7-11) 

The momentum fractions and the transverse momenta k±j^i are defined in Eqs. ()5.7|) 
and ()5.8j) . respectively, with (^j^ in Eq- The hatted momenta that also appear in the 

definition z-. ^ and are defined in Eq. ()7.13j) . 

Momentum mapping and phase space factorization 

The m momenta {pjm*'^'^*'' = {pi, . . . . . . ,pm+2} appearing in the matrix elements 

on the right hand side of Eq. ()7.5p are defined in two successive collinear mappings through 
an intermediate set of m + 1 momenta {p}m+i = {Pi, • • • ,Pkt, ■ ■ ■ ,Pm+2}, 

M ^ {]3}Si ^ {pyiF'''''^ , (7.12) 

or explicitly 

PL=r^(Pk+Pt-(^ktQn, P'n=r^Pn^ ^^k,t (7.13) 

1 - «fct 1 - ttfct 

and 

pTt, = T^^^P'kt+P'r-<ykt,Qn, P'n = TZ^Pn^ n^Tt,r, (7.14) 

^ ^ktr ^ ^ktr 

where a^t and are given by Eq. (j5.12p . the latter being defined on the hatted momenta. 

This iterated momentum mapping leads to exact phase space factorization in the iter- 
ated form 

d0,n+2(M;Q) = <lKmS''^''^-.QWiL'\vr,PTti^Q)]^^^^^ • (7.15) 

The one-parton factorized phase spaces [dp^^^^i{pk,Pkt]Q)\ and ['^Pi^m\pr^Pktf'^Qy\ 
given explicitly by Eq. (|5.16p . The graphical representation of this iterated momentum 
mapping and the corresponding factorization of the phase space is shown in Fig.lTUl 
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m+2 

Figure 10: Graphical representation of the iterated colhnear momentum mapping and the 
imphed factorization of the phase space. 

7.1.2 Iterated collinear— double collinear counterterm 
Counterterm 

Corresponding to the collinear limit of the double collinear subtraction is the countert- 
erm 

C.td;'2(M) = {Ma^^^^f — - (7.16) 

Skt S-f 

X (A^^^nM£^''=*))|PS,(^M, ^t,k, k^,k,u e)Pg{z.,, z,^, k^.,- e)\Mt\{p}t''^)) • 

The variables of the Altarelli-Parisi kernels, the momentum fractions and transverse mo- 
menta, are given by Eqs. ()5.7j) and ()5.8|) while the kernels themselves are recorded in 
Eqs. (Q-dSISI). 

Momentum mapping and phase space factorization 

The m momenta {p}™'^''^*^ = {pi, . . . ,pi~, . . . ,pkt, ■ ■ ■ ,Pm+2} in the matrix element on 
the right hand side of Eq. (j7.16j) are again given by an iteration of the collinear momentum 
mapping of Eq. (jS.lip , 

{P} ^ {P}^^li ^ {P}t''^ ■ (7.17) 
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Figure 11: Graphical representation of the iterated colhnear momentum mapping and the 
imphed factorization of the phase space. 



The new momenta are defined as 

P^, = r^(P'+Pr-^irQn, Pi^ = -^p^ n^lf, (7.18) 

1 - "if 1 - "if 

where the hatted momenta that also appear in the definition 2;- 2- • and fc_L ^'i- ^'^■1^1^ 
are defined in Eq. ()7.1H|1 . Again a^^, is given by Eq. ()5.12j) . This above iterated momentum 
mapping leads to exact phase space factorization in the following form: 

d0^+2(M; Q) = dM{p}t''^; g)[dpg(p.,Pj,; Q)W^^^M,Pku Q)] , (7.19) 

where the factorized phase space measures [dp['^^_^i{pk.,Pkt',Q)] and [dpi^^{pr,Pif,Q)] are 
given in Eq. (j5.16|) . The graphical representation of this iterated momentum mapping and 
the corresponding factorization of the phase space is shown in Fig.lTTl 



7.1.3 Iterated collinear— soft-collinear-type counterterms 



Counterterms 

The following three counterterms all use the same momentum mapping, that of the 
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coUinear-soft-collinear term Ci^^CS^^^.'^\{p}) and we present these together. We define 



Skt 

C,Ar-MCS^'t'?i{p}) = {^na.l^''f--^Tl (7.21) 



Skt 

= {Sna^n'n'^'f^Tl (7.22) 

ktf f,kt 

y<—{M^^\{p}^'''^)\P}lU^k,uz,^^^^^^^ 

Skt 

As usual Eq. ()5.7|) defines the momentum fractions and k±^k,t is defined in Eq. flfi.27|l . 
Momentum mapping and phase space factorization 

The set of m momenta {p}m^^^ = {pi, ■ ■ ■ ,Pkt, ■ ■ ■ ,Pm+2} {Pr is absent) that enter the 
matrix elements on the right hand sides of Eqs. ()7.20j) - (j7.22j) are constructed by applying 
the collinear momentum mapping of Eq. ()5.14|) followed by the soft mapping of Eq. ()5.26|) 
to the original set of m + 2 momenta {p} , 

{p} ^ ^ {p}Z'''^ , (7.23) 

which we have already discussed under Eq. (|6.29p and shown in Fig.lHl 



7.1.4 Iterated collinear— double-soft-type counterterms 
Counterterms 

The subtraction terms presented below use the same momentum mapping, thus they 
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are discussed together. We set 



1 

Skt 



)\\ (7.24) 



2e\2 



SktS 



-TICa 



X 



kt f 

^ r,fct / t Zth 



\ (1 -e)zk,tZt,k 



^fk 



-L.k,t 



' kt ,f 



Zt,k Zk,t 



'^±,k,t^ktf 



{0)fS^\{kt,kt)\\'2 



(7.25) 



2e\2 



kti 



X 



, kt 



^fk^^k,t 



+ Zk,tZt,k ,2 ^ 
^kt,f ^L,k,t^ktfJ 



\Mt\{p] 



n;(^fc«.fc«))|2 



(7.26) 



The momentum fractions and transverse momentum in Eqs. ()7.24|l - ()7.2(j|l are defined as 
usual by Eqs. ()5.7|) and ()5.8|) . In Eq. ()7.24p we use the notation 



^ kt kti 



(7.27) 



The hatted momenta are given by Eq. ()7.13p . 
Momentum mapping and phase space factorization 

The m momenta {pjm*'''"*^ = {pi, ■ ■ ■ ,Pm+2} (momenta andpj are absent) that appear 
in the matrix elements on the right hand sides of Eqs. ()7.24|) - ()7.26|) are obtained by applying 
the soft momentum mapping of Eq. ()5.26|) to the hatted set of momenta of Eq. ()7.13|1 . 



{p}^{pt!l^{pYl''''^ 



(7.28) 
(7.29) 



We have 

P''n = K[QAQ-Pkt)IHtWnlHt)- 

The hatted momenta are also used to define the momentum fractions z-. ^ in Eqs. (j7.24j) 
(|7.26p and are given in Eq. (j7.13p . 



The phase space factorization inherits the 'product' structure of the momentum map- 
ping and we get 



^+2(M;Q) 



Akt) 



{kt) 



(7.30) 
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Figure 12: Graphical representation of the coUinear momentum mapping followed by a soft 
one and the implied factorization of the phase space. 

The one-parton factorized phase space measures [(^Pi-^rn+iiPkiPkt] Q)] and [dp]^^ [p^t] Q)] are 
those in Eqs. (j5.16|) and ()5.28p . The graphical representation of this iterated momentum 
mapping and the corresponding factorization of the phase space is shown in Fig.fT2l 

7.2 Iterated soft counterterms 

7.2.1 Iterated soft— triple-collinear-type counterterms 

Counterterms 

The first three subtraction terms on the right hand side of Eq. ()7.3p are defined using 
the same momentum mapping, therefore these terms are presented together. They read 

^-{M^^\{p}t'^)\Pi%izi,,^^^^^ (7.31) 
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S{ir)t Zt .IT 

^^{Mt\{P]tHPf'fM.,,z,^^,k^^^^^^^ (7.33) 

The momentum fractions Zjj and Zj^^i are given respectively by Eqs. (j5.7p and ()6.5p while 
the transverse momentum k^i^is defined in Eq. (15. 8|) . The soft functions P^p^ j appearing 
in Eq. (j7.3ip were introduced in Ref. jl7j. As discussed below Eq. (j6.25|) . whenever j or / 
equals (ir) the eikonal factor appearing in Eq. ()7.32p evaluates to the expression given in 
Eq. (IfT^ . 

Momentum mapping and phase space factorization 

The m momenta {p}m'^^ = {pi, . . . ,pir, ■ ■ ■ ,pm+2} {pt is absent) entering the matrix 
elements on the right hand sides of Eqs. ()7.3ip - ()7.33|) are defined by the successive soft and 
collinear mappings, 

{P} ^ {pjJ^Vi ^ {P}t'^ > (7-34) 

where 

P':r = r^(Pi+Pr-^irQn, P'n = ^^P'n . ri ^ l^T (7.35) 

and the hatted momenta, that also appear in the definitions of z^-, z^- and k^--- in 
Eqs. (UnU-dlSni), are 

P'n = KIQ. {Q-Pt)/Xt]{p:/Xt) , n^t. (7.36) 

The expressions for aj- in Eq. (j7.35|) and At in Eq. (j7.36|) are given in Eqs. (j5.12|) and (j5.23|) 
respectively. This momentum mapping leads to an exact phase space factorization of the 
following form 

d<f>m+2{{p};Q) = d<jymmt^^ Q)[dpfZiPr^P^r■, Q)][dpfl^,{pu Q)] ■ (7.37) 

The collinear and soft one-parton factorized phase space measures [dp^'^^J^{pr,pir',Q)] and 

[dPi^m+iiPt'i Q)] are given respectively in Eqs. (j5.16p and (j5.28p . The graphical representa- 
tion of this iterated momentum mapping and the corresponding factorization of the phase 
space is shown in Fig.lT^ 

7.2.2 Iterated soft— double-soft-type counterterms 
Counterterms 
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Figure 13: Graphical representation of the soft momentum mapping followed by a coUinear 
one and the implied factorization of the phase space. 



The remaining four iterated soft counterterms in Eq. ()7.3|) are all defined using the same 
momentum mapping so we discuss them together. We have 



2e\2 



Oa 



1 Z: 



r,it 



1 '^i,rt 



2 Zi^t 2 z~. 



^it Zt.i S^,f. Z- ■ 



Sir ^r,i ^ J 2 



m 1 1 



(7.38) 
(7.39) 



SUfS„,,sT{{v}) = i8na.^^''rT^A'-^^^-^T^^\M^:^\{Pm\^ (7.40) 



2e\2 



I E '5^^(^)'5.'(t)|AlS),,,)(^.,)({p}(:'*))P (7.41) 



i,j,k,l 



i,k 



No new notation (besides the set of momenta to be defined below) has been intro- 

duced in Eqs. fl7.38|) - ()7.41|) . thus we limit the discussion of these terms to two comments. 
Firstly, note that in the equations above both r and t are gluons. Secondly, we emphasize 
that the discussion below Eq. (|6.25|) also applies to Eq. (|7.39|) . 
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Figure 14: Graphical representation of the iterated soft momentum mapping and the im- 
phed factorization of the phase space. 

Momentum mapping and phase space factorization 

The m momenta {p}m^^ = {pi, ■ ■ ■ ,pm+2] (both r and t are omitted) that appear in the 
matrix elements on the right hand sides of Eqs. ()7.38|) - ()7.41|) are obtained by two successive 
soft-type mappings (defined in Eq. ()5.22|l ). 

{P} ^ {P}^ti ^ {P}t;'^ > (7-42) 

where 

= K[Q, (g - pr)/\mj^r) , n^r. (7.43) 

The hatted momenta are given in Eq. (j7.36j) . Naturally, the iterative structure of the 
momentum mapping is inherited by the phase space factorization 

d0^+2(M; g) = d0^({p}(:'*); g)[dri^(pn g)][dM*L+i(pt; g)] , (7.44) 

where the soft one-parton factorized phase space measures [(^PilniPr'-, Q)] and [dpf^^j^-^pt] Q)] 
are both given by Eq. ()5.28p . The graphical representation of this iterated momentum 
mapping and the corresponding factorization of the phase space is shown in Fig. EI 

7.3 Iterated soft-collinear counterterms 

7.3.1 Iterated soft-collinear— triple-collinear-type counterterms 

Counterterms 
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The first two terms on the right hand side of Eq. ()7.4|l turn out to be defined using the 
same momentum mapping after a trivial reindexing. Let us define 

Sit Zt,i 

y<-{M^^\{p}t'^)\Pi^^^^^^^^^ (7.45) 

^if 

Skt Zt,k 

^^{M^^\{p}i''^)\Pili^lf,z,^^^ (7.46) 

Notice how a different indexing of CifSf-C^^f\{p}) is given above as compared to that 
appearing in Eq. ()7.4p . As already stated, this is convenient, because with this indexing 

exactly the same set of tilded momenta, {p}m'^\ appear in the matrix elements on the 
right hand sides of both Eqs. ()7.45p and ()7.46|) . The momentum fractions and transverse 
momenta are defined as usual via Eqs. (|5.7p and (|5.8|) . 

Momentum mapping and phase space factorization 

Exactly the same momentum mapping is used as for the iterated soft-triple-collinear- 
type counterterms. Sect. 17.2. 11 



7.3.2 Iterated soft-collinear— double-soft-type counterterms 
Counterterms 

The remaining five terms on the right hand side of Eq. ()7.4p are again defined using the 
same momentum mapping after reindexing some terms. We have 

C,ACS.,;tSir\{p}) = (8vra./.^^)^-^n-^T^|>l(;^)({p}(:'*))P, (7.47) 

^if ^f,i ^kt Zt,k 

C,ACkrtST\{p]) = {87ra.^^'r^'-^Tl^'^Tl\M^:^^m^'^)\\ (7.48) 

^kf ^f,k ^kt Zt,k 

C^ACkrt4r\{p}) = (87ra./.^^)^-^TL-^CA|A^(::)(W(:'*))r, (7.49) 

^kf ^f,k ^rt Zt,r 

CktStST\{p}) = -(Svras/i^^)^ ^ ^ ^5.Kr)f ^T^JA^;^),. ,({p}(;'*))P , 



. , 2 ^ SktZt,k 



(7.50) 
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(7.51) 

Here C^^SfCi^^fSl^''^\{p}) and C^j5(iS^j'°''({p}) are presented with a different indexing then 
in Eq. ()7.4j) so that all matrix elements on the right hand sides of Eqs. ()7.47j) - ()7.5H) appear 
with the same set of tilded momenta, {p}m^^- The momentum fractions have been defined 
in Eq. (j5.7|) . the eikonal factor in Eq. (j5.2ip . 

Momentum mapping and phase space factorization 

The momentum mapping used is identical with the mapping defined for the iterated 
soft-double-soft-type conterterms. This mapping and the corresponding phase space fac- 
torization is presented in Sect. 17.2.^ 



8 Cancellation of kinematical singularities 

In Sect.El we have shown that the subtraction terms collected in Eq. ()5.H) correctly regu- 
larize the kinematical singularities in the squared matrix element in the singly-unresolved 
regions of the phase space. Similarly, the subtraction terms collected in Eq. ()6.H) correctly 
regularize the kinematical singularities in the doubly-unresolved regions of the phase space. 
The purpose of the iterated counterterms is two-fold. These should cancel the kinemat- 
ical singularities of the singly-unresolved counterterms in the doubly-unresolved regions 
of the phase space and conversely, they should cancel the kinematical singularities of the 
doubly-unresolved counterterms in the singly-unresolved regions of the phase space. The 
structure of Eq. ()7.H) follows that of the candidate subtraction term Ai2|A1^^2p found in 
Ref . |17j , where it was shown that at the level of the factorization formulae the combination 
(A2 + Ai — Ai2)|A^||^l,_2p indeed regularizes the squared matrix element in all relevant un- 
resolved regions of the phase space. We show the structure of the cancellations graphically 
in Fig.[T31 

The first picture corresponds to the squared matrix element of the 771 + 2 final-state 
partons, while the following terms in the squared brackets correspond to the terms that 
build A2|7W^+2p, Ai|A1^^2p Ai2|A1^+2p) respectively. The factorized one- and two- 
parton factors correspond to Altarelli-Parisi kernels, with azimuthal correlations included, 
or eikonal factors, with colour-correlations included, or the combinations of these. In the 
first bracket we find all those terms that regularize the kinematical singularities of the 
squared matrix element in the doubly-unresolved regions of the phase space. In these 
phase space regions the terms in the second bracket, corresponding to Ai|A^^l,_2p, are 
regularized by terms in the third bracket, corresponding to Ai2|A^^!,_2p- In particular, if 
partons r and s become unresolved, than the two terms shown in the third line regularize 
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Figure 15: Graphical representation of the squared matrix element and its factorization 
formulae in the singly- and doubly-unresolved (soft and/or coUinear) limits. 

those two in the second line. 

In the singly-unresolved regions, the second bracket contains all terms, necessary to 
regularize the squared matrix element. Note however, that these terms also contain spurious 
singularities. For instance, when r is unresolved, then the first term in the second bracket 
regularizes the squared matrix element, while the second becomes singular. This term will 
be regularized by the second term in the third bracket. At the same time, the terms in the 
first bracket will lead to a strongly-ordered factorization formula, that will be regularized by 
the first term in the third bracket. If s becomes unresolved, then the role of the two terms 
in the third bracket interchanges: the first term will regularize the spurious singularity in 
the first term of the second bracket, while the second will regularize the term in the first 
bracket. 

In defining the full subtraction {yA.2 + A-i — Ai2)|A^m+2p5 "^^ keep the structure of 
(A2 + Ai — Ai2)|A^^!,_2p and replace the momenta in the squared matrix elements of m 
or m -|- 1 final-state partons in Fig.^J with tilded momenta, defined by different types of 
momentum mappings for the various terms. In order that the cancellations described in 
the previous two paragraphs take place it is crucial that these momentum mappings obey 
the following three conditions: 

1. The mappings used for defining the doubly-unresolved subtraction terms should be 
such that in the singly-unresolved regions of the phase space, the mapped momenta 
tend to the same limit as those of the iterated mappings, used for defining the terms 
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2. In those regions of the phase space, where only momentum becomes unresolved, the 
singly-collinear and soft momentum mappings, used in the definition of those singly- 
unresolved subtraction terms, in which momentum is factorized, tend to the same 
limit as the iterated mappings used in the definition of those terms in Ai2|A1m+2p; 
in which the first mapping factorizes and the second p^. 

3. In the doubly-unresolved regions of the phase space, the singly-collinear and soft 
momentum mappings, used in the definition of the ^i|A1m+2p term, tend to the 
same limit as the iterated mappings used in the definition of ^i2|Al||^+2p- 



The first of these conditions is necessary in order that the required cancellations be- 
tween v4.2|A^[||^2p ^i2|-^m+2p take place. The second requirement ensures that the 
cancellations between the first terms in the second and third bracket of Fig.^J happens 
when momentum becomes unresolved. Finally, the third condition is needed for can- 
celling all kinematical singularities of the v4.i|Al[||^2p terms by the corresponding terms in 
^i2|-A^m+2p i'^ the doubly-unresolved regions of the phase space. It is not difficult to check 
that the above requirements are fulfilled for all mappings. In particular, the third condi- 
tion follows from the construction of the iterative mappings, namely these are successive 
applications of the singly-collinear and/or soft mappings. Here we consider two illustrative 
examples of the first two conditions. 

In the case of the first condition, the least trivial is that the singly-collinear limit of the 
triply-collinear mapping is the same as that of two successive collinear mappings, defined 
by Eqs. (j7.12j) - (j7.14j) . as shown graphically in Fig.^l In the limit when momenta p^ and 
pf are collinear, using Eqs. (|6.5p and we find 

Zk,ts ^ ^k^{kt),s ! Zt,ks ^ ZtZ{jtt),s , Zs,kt ^ (8-1) 

PkWpt PkWpt PkWpt 

and 

^±,fc,ts ^kk'l (kt),s ' f^±,t,ks ^tk'l,(kt),s ' k±,s,kt ^±,s,(kt) ' (8-2) 

PkWpt PkWpt PkWpt 

i.e., both the momentum fractions and the transverse momenta tend to the limit of the 
corresponding variables of the iterative collinear mapping. As a result, the kinematics 
defined by the two mappings tends to the same limit and the cancellation of the singularities 
takes place. 

For the second requirement, let us consider for example, the soft limit — 0. The 
singly- unresolved counterterms depend on p^, while the iterated terms depend on p^. The 
singly-unresolved mappings, that lead to these different momenta in the two cases, are 
linear in p^, therefore, — > in the soft limit, when p^ — *■ 0. It follows that ^ 1 when 
p^ — > so in the soft limit the hatted momenta tend to those with tilde. In particular, 
Ps ~* Ps if Ps 0' so the kinematics become identical and the cancellation takes place. 
The graphical representations of the soft limits of the two mappings are shown in Fig. 1171 
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m+2 m+2 m+2 

Figure 16: Graphical representation of the colhnear hmit of the triply-colhnear and and 
iterated singly-colhnear mappings, (kt) means the momentum p'j^ + in the colhnear 
direction in the colhnear limit. 



11 i 




m+2 m+2 m+2 



Figure 17: Graphical representation of the soft limit of the singly-unresolved and iterated 
momentum mappings. 
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9 Numerical results 



In Sections El IHl and [3 we have defined explicitly all subtraction terms in Eqs. ()4.H|1 - ()4.5|1 

that are necessary to make dcx^^j integrable in ci = 4 dimensions. In order to further 
demonstrate that the subtraction terms indeed regularize the cross section for doubly-real 
emission, we consider the non-trivial example of the contribution of the e+e~ —>■ qqggg 
subprocess to the moments of three-jet event-shape variables thrust (T) and C-parameter, 
when the jet function is a functional 

J„(pi, . . . O) = 6{0 - Osipu ■ . .,Pn)) , (9.1) 

with 0^{pi, . . . ,Pn) being the value of either r = 1 — T or C for a given event (pi, . . . ,Pn)- 
The e~^e~ — qqggg subprocess gives rise to the most general kinds of kinematical singu- 
larities and colour structures. The only additional complication present in the four-quark 
subprocess e~^e~ qqQQg is the identical flavour contribution that does not require any 
addition to the subtraction scheme. 

Starting from randomly chosen phase space points and approaching the various singly- 
and doubly-unresolved regions of the phase space in successive steps, we have checked 
numerically that the sum of the subtraction terms has the same limits (up to integrable 
square-root singularities) as the squared matrix element itself. 

The perturbative expansion to the n^^ moment of a three-jet observable at a fixed scale 
fi = Q and NNLO accuracy can be parametrized as 

(O") = /"dOO"— ^ = 
J ctq dO 

^) -^o ' + (^) ' «o ' + (^) ' • f^-^) 
where according to Eq. ()3.4j) . the NNLO correction is a sum of three contributions 

= 41 + C^ol + C^ol ■ (9.3) 

Carrying out the phase space integrations in Eq. (j3.5p . we computed the five-parton contri- 
bution Cq}^{0) as defined in this article. The predictions for the first three moments of r 
and the C-parameter, obtained using about ten million Monte Carlo events, are presented 
in Tabled These numbers are unphysical, and given only to demonstrate that the (m + 2)- 
parton NNLO cross section defined in this paper is finite. In particular, the relatively small 
negative values simply indicate that the subtractions altogether subtract slightly more then 
the full doubly-real cross section. If needed, the colour decomposition is straightforward. 
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Table 1: The moments C^l and C^ 



n 






1 


-(9.27 ±0.33) 


10 


-(3.44 ±0.13) 


10^ 


2 


-(0.31 ±0.04) 


10 


-(1.41 ±0.03) 


10^ 


3 


-(0.20 ±0.01) 


10 


-(0.63 ±0.18) 


■ 10 



10 Conclusions 



In this paper we have set up a subtraction scheme for computing NNLO corrections to QCD 
jet cross sections to processes without coloured partons in the initial state. The scheme is 
completely general in the sense that any number of massless coloured final-state partons 
(massive vector bosons are assumed to decay into massless fermions) are allowed provided 
the necessary squared matrix elements are known. It is also general in the sense that it is 
algorithmic in a straightforward manner, therefore, the generalization to N^LO accuracy 
does not require new concepts. Each step of the computation can in principle be extended 
to any order in perturbation theory, which is useful in setting up parton shower algorithms 
that can be matched to fixed-order approximations naturally. 

Three types of corrections contribute to the NNLO corrections: the doubly-real, the 
real-virtual and the doubly-virtual ones. Here we have constructed the subtraction terms 
for the doubly-real emissions; those to the real-virtual corrections will be presented in a 
companion paper. By rendering these two contributions finite in d = 4 dimensions, the 
KLN theorem ensures that for infrared safe observables adding the subtractions above to 
the doubly-virtual correction, that becomes also finite m d = A dimensions. 

The subtraction terms for the doubly-real corrections presented here are local m d = 
A — 2e dimensions and include complete colour and azimuthal correlations. The expressions 
were derived by extending the various singly- and doubly-unresolved limits of QCD squared 
matrix elements over the whole phase space, which was achieved by introducing momentum 
mappings which define exactly factorized phase-space measures. Although the number of 
subtraction terms is rather large, the implementation of the scheme is not so complicated, 
because only five different types of phase-space mappings have to be defined, all other 
mappings being obtained by employing those five basic ones iteratively. 

In order to demonstrate that the subtracted cross section is indeed integrable, we have 
computed the corresponding contributions to the first three moments of two three-jet event- 
shape observables, the thrust and the C-parameter. 
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